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We report the first detection of the lowest CO transition in a sub-millimetre 
bright galaxy and extremely red object (ERO) at z — 1.44 using the Very Large 
Array 1 . The total J = 1 - line luminosity of ERO J164502+4626.4 is (7±1) x 
10 10 Kkms" 1 pc 2 , which yields a total molecular gas mass of ~ 6 x 10 10 M & . We 
also present a map of the 850-/iin continuum emission obtained using SCUBA, 
from which we infer a far-IR luminosity and dust mass of Lfir ~ 9 x 10 12 L® and 
M d ~ 9 x 10 8 M Q . We find tentative evidence that the CO and sub-mm dust 
emission is extended over several tens of kpc. If confirmed by high-resolution 
imaging, this implies that ERO J164502+4626.4 is not simply a high redshift 
counterpart of a typical Ultra Luminous Infrared Galaxy (ULIRG). 

Subject headings: galaxies: structure - galaxies: individual: (ERO J164502+4626.4) 
- galaxies: ISM 



1 The Very Large Array (VLA) is operated by the National Radio Observatory, which is a facility of the 
National Science Foundation, operated under cooperative agreement by Associated Universities, Inc. 
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1. Introduction 

Observations of CO provide one of the most powerful methods of probing the interstellar 
medium (ISM) in galaxies - i.e. determining the amount of molecular gas available to fuel 
star formation and accretion onto active galactic nuclei (AGN), the two processes believed 
to generate the large far-infrared (far-IR) luminosities of star-forming galaxies (Sanders et 
al. 1991). The first detection of CO at high redshift (Brown and vanden Bout 1991) revealed 
its potential to trace metal-enriched molecular gas in the early Universe. CO is thus key to 
our understanding of galaxy formation and evolution. 

Extremely red objects (EROs - usually defined as galaxies with R—K > 5.3, e.g. Moriondo, 
Cimatti & Daddi 2000) constitute a bi-modal population - a combination of dusty, starburst- 
ing systems (Dey et al. 1999) and evolved ellipticals (Dunlop et al. 1996). While both types 
of objects appear similar in the optical/near-IR, the detection of sub- millimetre continuum 
emission and/or CO unambiguously pinpoints the dusty EROs. A significant fraction of 
the sub-mm-selected population of high-redshift dust-enshrouded starburst galaxies is be- 
lieved to be associated with EROs (Smail et al. 1999). Since the advent of large- format 
bolometer arrays such as SCUBA (Holland et al. 1999) and MAMBO (Kreysa et al. 1998), 
over a hundred sources have been detected in (sub) millimetre surveys (Smail et al. 1997; 
Hughes et al. 1998; Bertoldi et al. 2001). In cases where optical/IR counterparts of sub-mm 
galaxies are available, it is found that the sources typically are distorted, multi-component 
systems with one or more components being an ERO (Ivison et al. 2002). The sub-mm 
population has a median redshift of z ~ 2.4 (Chapman et al. 2003), and is widely believed 
to comprise the progenitors of present-day spheroids and massive ellipticals. At present, 
however, the molecular gas content of sub-mm selected galaxies remains largely unknown, 
and a systematic inventory of the molecular gas and its properties in these sources is required 
in order to properly address the question of their typical mass and evolutionary status. Of 
particular importance is the ratio of gas and dynamical masses in such objects, since it can 
help determine whether sub-mm galaxies are merely high-redshift analogues of local ultra 
luminous infra-red galaxies (ULIRGs) or massive, large-scale galaxy-formation events. The 
difficulties arise mainly because sub-mm selected galaxies have proved extremely difficult to 
identify in the optical (cf. Ivison et al. 1998, 2000). As a result only a handful of sub-mm 
selected galaxies have been detected in CO to date (Frayer et al. 1998, 1999; Ivison et al. 
2001; Downes & Solomon 2003; Genzel et al. 2003; Neri et al. 2003) 

ERO J164502+4626.4 (Graham & Dey 1996) was amongst the first EROs discovered 
by Hu & Ridgway (1994). Deep near-IR and optical spectroscopy of J164502 put it at a 
redshift of z — 1.44 (Graham & Dey 1996; Dey et al. 1999), and the detection of the Ha 
line and the [O II]A3726,3729 doublet in its spectrum suggested that J164502 is an actively 
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star-forming galaxy, possibly containing an AGN, and not an evolved elliptical at z = 2 — 4 
as initially suggested by Hu & Ridgway (1994). The presence of large amounts of dust was 
inferred from sub-mm observations by Cimatti et al. (1998) and Dey et al. (1999), the latter 
detecting J164502 in the continuum at 450, 850 and 1350 fim. These observations together 
with detections of CO J = 2 — 1 and J = 5 — 4 by Andreani et al. (2000), which revealed the 
presence of large quantities of molecular gas, unamibiguously demonstrated that J 164502 is 
a gas-rich dust-enshrouded starburst galaxy. 

Here we present observations of the CO J = 1 — emission from J 164502 using the 
VLA, providing the first detection of the lowest, and thus least excitation biased, CO line 
in a galaxy typical of the sub-mm population. We also present the 850-/zm SCUBA map 
of this source. Throughout this paper we have assumed H = 65kms _1 Mpc" 1 , Vt M = 0.3 
and = 0.7. In this cosmology the luminosity distance of J164502 is 11.2 Gpc and 1" 
corresponds to 9.1 kpc. 

2. Observations and data reduction 

2.1. Very Large Array observations 

At z = 1.44 the CO J = 1 - line from J164502 falls within the VLA's Q-band at 
0.7 cm. The line widths of the CO J = 5 - 4 and J = 2-1 lines were ~ 400 km s -1 FWHM 
(Andreani et al. 2000), and a similar line width could be expected for the J = 1 — line. 
In order to cover the entire CO line, two adjacent 50 MHz IF pairs (left and right circular 
polarisation) were centred at 47.235 and 47.285 GHz giving 635 km s -1 of coverage. Mea- 
surements of J164502 in the radio have yielded S^ocm < 300 /iJy and S3. % cm = (35 ±11) A*Jy 
(Frayer 1996). Fitting a power law spectrum to these data points yields a spectral index 
of a > —1.25. Adopting a more conservative spectrum with a = —0.7, we find that the 
extrapolated synchroton flux at 0.7 cm amounts to no more than 11 /iJy. Furthermore, the 
dust contributes by less than 10 /iJy at these wavelengths (Dey et al. 1999). Hence, any 
contribution from continuum emission is neglible and no continuum subtraction was nec- 
essary. Observations were obtained in D and C configurations (see Table 1). Due to bad 
weather, all the C-configuration data were of poor quality and discarded, leaving a total of 
11.6 hr of observing time on source. The primary calibrator, 3C286, was observed at the 
beginning and end of each observing run in order to fix the absolute flux density scale. Every 
3-4 minutes the antennas were pointed toward a phase calibrator to ensure phase coherence 
throughout the run. This fast-switching technique resulted in a residual phase rms, after 
calibration, of less than ~ 20°. Calibration and data reduction was done using standard 
recipes in the NRAO AIVS Cookbook. D-configuration data taken at different times were 
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combined using dbcon, weighting each dataset by its total gridded weight. 

The astrometrical uncertainty in the final CO map is given by a 2 = a 2 al + a 2 es + cr 2 SNRl 
where a ca i is the uncertainty of the assumed position of the phase calibrator, a res is the 
residual rms phase error after the phase correction, and cjsnr is the positional error due to 
the presence of thermal noise. In our case the latter is by far the most dominant source of 
error and is given by a SNR ~ (FWHM/2) / SNR. Adopting FWHM ~ 3" and SNR ~ 7 
(see Fig. 1) we find the total positional error in the final CO map to be ~ 0.2". 



2.2. Submillimetre observations 

Imaging data at 450 and 850-/zm were obtained for J164502 using the Submillimetre 
Common User Bolometer Array (SCUBA - Holland et al. (1999)) during periods of ex- 
ceptionally low opacity at the James Clerk Maxwell Telescope (JCMT 2 ) on Mauna Kea, 
Hawaii. 

Data for J164502 were acquired during 1998 July 28-30 and 1998 September 03-05, a 
total integration time of 28.8 ks using a 45" east-west nod/chop. In brief, data were flat- 
fielded, corrected for atmospheric losses, edited both automatically and by hand, plotted 
onto an astrometric grid and flux calibrated using observations of planets and secondary 
calibration sources. The map was convolved with a Gaussian resulting in a resolution of 14". 



2.3. Near-IR Observations 

Near-IR K-band imaging was obtained of J164502 in June 2002 using the United King- 
dom IR Telescope's (UKIRT 3 ) UFTI imager, a 1024 2 HgTeCd array with 0.091" pixels. Nine 
individual frames taken during good to moderate seeing conditions were co-registered and 
combined yielding a total exposure time of 2700s. The seeing in the final image is ~ 0.8". 

Using 5 stars also detected in the USNO A2.0 catalogue (Monet et al. 1998), the UKIRT 
image was calibrated on to the USNO reference frame to within 0.2" rms. The CO map was 



2 The JCMT is operated by the Joint Astronomy Centre on behalf of the United Kingdom Particle Physics 
and Astronomy Research Council (PPARC), the Netherlands Organisation for Scientific Research, and the 
National Research Council of Canada. 

3 UKIRT is operated by the Joint Astronomy Centre on behalf of the United Kingdom Particle Physics 
and Astronomy Research Council (PPARC) 
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aligned to the K-band frame using the tasks GEO and hgeom in AZVS resulting in a relative 
positional accuracy between the two of ~ 0.3". 



In order to sensitively search for possible faint extended CO J = 1 — emission, at the 
expense of high spatial resolution, we used the compact D-configuration data tapered down 
using a Gaussian with a full-width of 60 kA at 30 per-cent amplitude, yielding a synthesized 
beam of 3.1" x 2.8". The resulting CO map overlaid on the UKIRT K-band image is shown in 
Fig. 1. CO J = 1 — emission is detected at the 7a significance level in J 164502. Assuming 
that we cover the entire J = 1 — line with the two IFs, and that we have rectangular 
passbands, we estimate the velocity-integrated flux-density using 



where Au IF is the width of the IF and Sco = (Si + 5 , 2)/2 is the average of the flux density 
measured in the two IFs. Using the task imean to get the flux density of the source in the 
two IF maps we find S C o — 1-0 ± 0.2 mJy which yields a velocity-integrated flux density of 
0.6 ± 0.1 Jykms -1 , where we have used v b s = 47.26 GHz and an IF bandwidth of 45 MHz 
instead of 50 MHz due to bandpass rollover. 

As first noted by Dey et al. (1999) using high-resolution HST WFPC2 images, J164502 
in the optical has a reflected S-shaped morphology with two bright knots at each end. A 
comparison with Keck K-band imaging revealed that the bulk of the near-IR emission comes 
from the region between the two bright knots located ~ 0.4" north of the brightest optical 
knot (Dey et al. 1999). As seen from Fig. 1, the peak intensity in the CO J = 1 — map 
coincides with the peak in the UKIRT K-band image both of which are coincident with 
the region of low optical emission found by Dey et al. (1999). A similar offset between 
the HST position and the CO J = 2 — 1 emission was reported by Andreani et al. (2000). 
While the near-IR emission traces the dust-enshrouded stars in J164502, the optical emission 
corresponds to regions of low extinction, and it is therefore not surprising to see such an 
offset between the two types of emission. Since gas should be a good tracer of dust one 
would expect the CO to coincide with the K-band emission, as it appears to do within the 
limits of our astrometric errors. 

A Gaussian fit to the CO data in the image plane yields a deconvolved source size of 
4.5" x 3.0", suggesting that even in the D-configuration we have resolved the CO emission, 
albeit marginally. Different weighting schemes and tapering functions were tried in order to 
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Results 
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see the effects on the extended emission. Using a less restrictive tapering function (FWHM= 
200kA) and a robust parameter of 2 instead of 5 (natural weighting), thereby increasing 
the weight of the long baselines, the source still appeared extended although at a lower 
significance due to the larger noise on the longest baselines (see Figure 1). In order to confirm 
or dismiss the reality of the extended emission we shifted the peak of the CO emission to 
the phase center, and then plotted the binned real and imaginary parts of the complex 
visibilities as a function of baseline distance, Fig. 1. The complex visibilities are the Fourier 
components of the source brightness distribution. Hence, for a Gaussian source brightness 
distribution the real part of the complex visibilities is a Gaussian while the imaginary part 
will be zero since a Gaussian is an even function. The dashed curves in Fig. 1 represents 
a least-squares fit of a circular Gaussian to the visibilities. The position of the Gaussian 
was fixed to the phase center during the fit and only the amplitude and width were allowed 
to vary. We find that the real part of the visibilities is consistent with a Gaussian source 
model with a FWHM = 3.9" in agreement with the Gaussian fit in the image plane. The 
imaginary part is consistent with zero at all baselines within the limits of our residual rms 
phase error of ~ 20°. The phase error would give rise to a dispersion around zero of 
a(Im(V)) = a(Ssin0) ~ Sa(<f>) =1.0 mJy (20°/180°)7r ~ 0.3 mJy which is close to the 
observed scatter. 

The 850-yitm SCUBA map of J164502 is shown in Fig. 2, and is seen to contain at least 
two statistically significant sources. A 7er emission feature showing in the map is detected at 
a position which coincides with the near-IR and CO-detections of J164502, and we take this 
to unambiguously be the 850-/xm SCUBA detection of J164502. Another source, marked as 
J164502-SMM1 in Fig. 2, detected at the 5a level is seen 50" to the south, with faint emission 
at the 4<T level also seen just ~ 20" south of J164502. We estimated the 850-/zm flux density 
in four different ways. Two estimates were made by measuring the flux within apertures of 
20" and 16" diameter. A third measurement was made by fitting a Gaussian to the emission, 
and a fourth by fitting a Gaussian with a fixed FWHM of 15" to the emission. The values 
derived from each of these measurements are given in Table 2 for J164502 as well as J164502- 
SMM1. For J164502 we find S 850 ^ m = 8 ± 2mJy where we have adopted the average of the 
four different measurements. While this is in good agreement with the 850-/zm measurement 
by Cimatti et al. (1998), it is larger than the flux quoted by Dey et al. (1999). Both of these 
measurements were obtained using SCUBA in its photometry mode but the latter value is 
based on twice as large a data-set as the former, taken in excellent weather conditions and 
is therefore the most reliable estimate. If this is the case, the large discrepancy between the 
flux density derived from our map and the value derived from photometry-measurement by 
Dey et al. (1999) could be due to extended emission from J164502 which would have been 
missed by the single pixel measurement. 
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From the 850-/xm map in Fig. 2 it is seen that the emission from J164502 does in fact 
appear to be extended, with emission to the NE. The reality of this extended emission is 
strengthened by the high- resolution Keck K-band image obtained by Dey et al. (1999) which 
shows an emission feature extending in the same direction, albeit on a smaller scale of > 1" 
(see Fig. 2 of Dey et al. (1999)). However, the fact that the stars traced by the K-band 
emission is extended on scales of >9 kpc means that cold dust is distributed on at least the 
same scales, and likely on much larger scales. Fitting a Gaussian to the SCUBA image yields 
a source size of 15.1" along the major axis, which is comparable to the 14" resolution of the 
map, though very tentatively suggestive of a slight extension. In Fig. 2 we compare the 
azimuthally averaged radial profiles of J164502 and J164502-SMM1 with the point spread 
function (PSF) of SCUBA at 850-/im derived from a beam map of a bright blazar obtained 
during the same observing run. The radial profile of J164502 agrees with the PSF out to 
r ~ 20" beyond which it displays excess emission over that of a point source. This bump 
in the radial profile is due to the 4a emission feature 20" south (~ 180 kpc) of J164502. 
The reality of this emission feature is questionable but if real its close vicinity to J 164502 
could mean that it is a system in the process of merging with J164502. While there is thus 
tentative evidence that J164502 is extended, the large beam size of SCUBA does not allow 
for a firm conclusion on this issue. 

For comparison we note that J164502-SMM1 is consistent with the PSF at all radii. 
The SCUBA 850-/im beam is too large to make a reliable identification of J164502-SMM1 
with an near-IR/optical source. If J164502-SMM1 is a dusty system at the same redshift as 
J164502 one might expect that it too would shine in CO. No CO emission was detected at 
the position of J164502-SMM1, however, and we set an upper limit on the velocity- integrated 
line flux using the equation Sco^v < 3a(5vdvY^ 2 , where a is the rms noise, Sv the velocity 
resolution, and dv the line width which we assume to be equal to the J = 2 — 1 line width, 
i.e. ~ 400 km s -1 (Andreani et al. 2000). The measured rms noise in the image is o ~ 
0.1 mJy beam -1 which yields Sco^v < 0.2 Jykms -1 for a velocity- coverage of 635 km s -1 . 
The failure to detect CO-emission from J164502-SMM1 does not rule out a association since 
cluster velocity dispersions can go up to > 1000 km s -1 , and in some cases even closely 
associated CO-emitting regions can have velocity differences of ~ 1000 km s -1 as was seen 
in the high redshift radio galaxy 4C 60.07 (Papadopoulos et al. 2000). The probability of 
finding by chance a source with a flux density brighter than or equal to 6 mJy within r = 50" 
from J164502 is given by P = 1 — e -7rr2jv , where N is the average surface density of sources 
with 5*850^111 > 6mJy (Downes et al. 1986). Adopting a surface density of N ~ 400 deg -2 
(Borys et al. 2003) one finds that P = 0.22. On that basis we conclude that J164502-SMM1 
is most likely unrelated to J164502. 
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4. Analysis & Discussion 

Table 2 summarises the CO and sub-mm observations of J164502 and J164502-SMM1 
as well as the physical quantities derived from them. 

4.1. Molecular gas and dust in J 164502 

For an observed velocity-integrated line flux density of the J + 1 — > J CO line, the 
intrinsic CO line luminosity is given by 

c 2 D 2 f 

Lc °™ = 2k^;rT^ y Av s — dv ' (2) 

where D L is the luminosity distance, and u J+ i^ is the rest-frame frequency of the CO 
transition. Inserting astrophysical units, the equation reads 

L< 00J+1J = 2.43 x 10.(J + !)-(! + , (3) 

where L'co 7+ i 7 is the so-called pseudo-luminosity which is measured in units of Kkms -1 pc 2 . 
The observed CO J = 1 — line flux density for J164502 implies an intrinsic CO luminosity 
of L' co (l - 0) = (7 ± 1) x 10 10 Kkms^pc 2 . From the CO J = 1 - line the molecular 
gas mass can be found using the scaling relation M(H 2 ) = X co L co (l — 0) where X co ~ 
5M (Kkms -1 pc 2 ) -1 is the standard Galactic CO-H 2 conversion factor. The standard 
conversion factor, X C o has been calibrated using Giant Molecular Clouds (GMCs) in our 
Galaxy (e.g. Strong et al. 1988). However, the conditions of the interstellar medium in 
active starbursting galaxies such as ULIRGs is markedly different from that in our Galaxy, 
and applying the Galactic conversion factor would result in seriously underestimated gas 
masses. In particular, the concept of isolated and virialised GMCs breaks down, and instead 
the bulk of the CO emission is found to originate from a warm diffuse phase. In effect, this 
yields a lower conversion factor of X co = 0.8 M (Kkms -1 pc 2 ) -1 , more appropriate for 
such extreme environments (Downes & Solomon 1998). It is reasonable to assume that in 
high-redshift starbursts similar conditions prevail and the aforementioned value of Xco is 
also adopted for J164502. In doing so we estimate the molecular gas mass present in J164502 
to be M(H 2 ) = (6 ± 1) x 10 10 M . This is in agreement with the J = 2 - 1 detection by 
Andreani et al. (2000), which yields a molecular gas mass of ~ 4 x 10 10 M when converted 
to the cosmology adopted in this paper. 

The CO J = 1 — line luminosity and the implied molecular gas mass of J 164502 is 
somewhat lower than those found in the quasars APM 08279+5255 and PSS 2322+1944, and 
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in the high redshift radio galaxy 4C60.07 (L' co (l - 0) ~ 10 11 Kkms" 1 pc 2 and M(H 2 ) ~ 
10 11 M Q - see Papadopoulos et al. 2001; Carilli et al. 2002; Greve et al. in preparation). 
This is not surprising since these objects are among the most luminous systems in the 
Universe, and might therefore be expected to contain more gas. Solomon et al. (1997) 
observed CO J = 1 — in a sample of 37 local ULIRGs out to a redshift of z — 0.3, and found 
an average J = 1 — luminosity of L' co (l — 0) ~ 8 x 10 9 Kkms -1 pc 2 which is roughly the 
luminosity of Arp220 and Mrk231. The scatter on this result was only 30 per cent. A more 
recent survey of CO J = 1 — observations of a complete sample of 60 ULIRGs selected from 
the SCUBA Local Universe Galaxy Survey (SLUGS - Dunne et al. (2000)) yielded similar 
results (Yao et al. 2003). Hence, we find that J164502 has a CO J = 1 — luminosity and 
consequently a molecular gas mass which is about an order of magnitude larger than is found 
for local ULIRGs. 

The dynamical mass of the system can be calculated from the observed size and line 
width of the source and assuming the gas is distributed in a disk with diameter L in Keplerian 
rotation, in which case it can be shown that the dynamical mass is given by 

M _ av fwhm l 1 (4) 

LOldLx sin % 

11/- , n 9 / AVpwHM \ 2 ( R \ I , , /_x 

= 1.16 x 10 M T x — o~Mq , (5) 

V100kms~V V k P c /sin 2 i °' v; 

where % is the inclination angle of the disk, and a d is a correction factor of order unity 
(Bryant & Scoville 1996). Adopting a line width of AVfwhm = 400 km s" 1 (Andreani et 
al. 2000), a source size of 9 < V4.5" x 3.0" = 3.7", which corresponds to a maximum linear 
diameter of L ~ 34kpc at a redshift of z — 1.44 we estimate the enclosed dynamical mass 
within the CO-emitting region to be M^ n < 6.3 x 10 n (sini)~ 2 M G . The inferred ratio of 
the molecular-to-dynamic mass for the system as a whole is then M(H 2 )/Mdy n > 0.1 sin 2 i. 
Hence, geometrical factors aside, the amount of molecular gas accounts for at least 10 per 
cent of the total dynamical mass within the CO-emitting region. The dynamical mass should 
be considered an upper limit on the total amount of molecular gas present in J164502, and 
the two mass estimates would only coincide if the CO emission was concentrated within the 
inner ~ 3 kpc. Due to the higher spatial resolution and the fact that the J = 1 — line 
provides an unbiased estimate of M(H 2 ), our constraint on the molecular-to-dynamical mass 
ratio should be an improvement over that of Andreani et al. (2000). A gas mass fraction of 
> 10 per cent is compatible with the range of gas mass fractions found for local ULIRGs 
(Downes & Solomon 1998) yet, unlike the latter, it is distributed over significantly larger 
scales. 



The extrapolated synchroton radiation flux at 850- jum is less than 3 /xJy, and it is there- 
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fore safe to assume that the observed 850-/xm flux from J164502 is dominated by thermal 
dust emission, with the radio synchroton emission contributing a negligible amount. The 
dust mass can therefore be estimated using 



where v rest = ^ & s (l + z ) is the rest-frame frequency, T cm b(z) is the CMB temperature at 
redshift z, and Kd oc u 13 is the dust absorption coefficient. The emissivity index, /3, depends 
on the dust temperature, T d , (Dunne et al. 2000) both of which are poorly constrained for 
high-redshift objects - typical values are f3 — 1 — 2. Dey et al. (1999) found the spectral 
energy distribution (SED) of J164502 to be well described by an optically thin modified black 
body law with a dust temperature of T d = 40 K and an emissivity law of f3 — 1.5, and we 
adopt those values here. In doing so we estimate the dust mass to be M dU st — 9 x 10 8 M Q , 
where we have used a dust absorption coefficient of Kd{v r est) — 0.11(z/ rest /353 GHz)^ m 2 kg -1 
after Hildebrand (1983). The uncertainty in k d is large, so we present this estimate for 
comparison only. We can estimate the total far-IR luminosity by integrating the thermal 
spectrum 



which yields a total far-IR luminosity of Lpm — 9 x 10 12 L Q . This almost puts J164502 in 
the class of hyperluminous infrared galaxies (HyLIRGs) which have Lpm > 10 13 L & . The 
dust masses and far-IR luminosities typically found in local ULIRGs are M dU st — 10 7 ~ 8 M Q 
and Lfir > 1O 12 L (Sanders & Mirabel 1996; Dunne et al. 2000), i.e. nearly an order of 
magnitude smaller than what we find for J164502. It is comparable, however, to what is 
found in high-redshift QSOs and HzRGs (Omont et al. 2003; Archibald et al. 2001) and, 
perhaps more importantly, similar to the dust masses and far-IR luminosities found in sub- 
mm-selected dust-enshrouded starbursts at high redshifts (Ivison et al. 1998, 2000). From 
the above we estimate the gas-to-dust mass ratio of J164502 to be M(H 2 )/Md us t — 67 which 
is well within the range of 50-100 found for local ULIRGs. Note, however, that significant 
uncertainty is attached to the normalisation of Kd, and values differing by as much as a 
factor of two has been reported (Draine & Lee 1984; Mathis & Whiffen 1989), which in turn 
could lead to a similar change in the dust mass but not in Lpm which is independent of the 
normalisation value of the dust absorption coefficient. Using (3 = 2.0 instead of 1.5 would 
decrease the dust mass by nearly a factor of two, but increase L FIR by a similar amount. 

If the bulk of the far-IR emission can be ascribed to starburst activity the corresponding 
star formation rate is given by 






(7) 



SFR ~ L 




>imf$sb L Q M q yr 



-i 



(8) 
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where 5jmf ~ 1—6 is a function of the initial mass function, and Ssb is the fraction of the FIR 
emission which is heated by the starburst (Omont et al. 2001). Assuming a conservative value 
of 

$imf — 1) we estimate the starformation rate in J164502 to be SFR ~ 9005sb M Q yr . 
It is possible that J 164502 harbours an AGN in its center, so a significant fraction of the far- 
IR luminosity could be due to dust being heated by the AGN and not by the starburst. The 
derived value for the SFR could therefore be overestimated. However, the narrow linewidths 
seen in optical and NIR spectra of J164502 (Dey et al. 1999) favours young hot stars over 
an AGN as the main source of energy. Furthermore, as pointed out by Dey et al. (1999), 
J164502 deviates from the 60 /iin - 6 cm relation obeyed by local star-forming galaxies by 
having almost an order of magnitude less radio emission at 6 cm, which is unlikely if an AGN 
dominated the energetics of the system (Dey & van Breugel 1994). In addition, it is found 
that the bulk of the far-IR emission from local ULIRGs is powered by starburst activity even 
though an AGN is present (Downes & Solomon 1998). Finally, if the sub- mm emission is 
extended it would strongly suggest that the main power source for the far-IR emission is a 
massive starburst since it is difficult to imagine the AGN heating the dust on scales of tens 
of kpc. The large quantity of molecular gas revealed by the CO detections in this system 
could provide the necessary fuel for such a massive starburst for a period of ~ 10 7 years. 

The efficiency of star formation should be measured relative to the amount of molecu- 
lar gas available to form stars. Such a measure of the star-formation efficiency is the rate 
of star formation per solar mass of molecular hydrogen, i.e. L FIR /M(R2)- For J164502 
we find L FIR /M(R 2 ) — 150 L M _1 . Probably a better gauge of the star-formation effi- 
ciency is the continuum-to-line ratio L FIR / L' co (l — 0) since it is independent of X co . Giant 
Molecular Cloud (GMCs) in our Galaxy typically have values of ~ 15 and similar ratios are 
found in nearby spirals (Mooney & Solomon 1988). Starburst galaxies and ULIRGs have 
Lpi R j L' co (l — 0) ratios which are 10 times higher than this, ranging in values from 80 to 
250 with a median of 160 L G ( Kkms" 1 pc 2 )" 1 (Solomon et al. 1997). For J164502 we find 
L FIR /L' co (l — 0) ~ 129 L G ( Kkms -1 pc 2 ) -1 , i.e. in line with what is found for local ULIRGs. 
Recent studies of large samples of ULIRGs at low and intermediate redshifts have shown 
that the L FIR /L' co (l — 0) ratio increases with increasing L FIR (Young et al. 1986; Tutui 
et al. 2000; Yao et al. 2003). Such behaviour can be explained if the more far-IR luminous 
galaxies, in addition to having more dust and gas which would just continue the linear re- 
lation between L F i R and L' co (l — 0), have higher dust temperatures due to extra heating 
by an AGN or a higher star-formation efficiency. Assuming that there is no contribution 
to L F j R from an AGN, the derived Lpi R / 'Z/ co (l — 0)-ratio should probably be taken as a 
lower limit on the star-formation efficiency since in the extreme starburst regions dominating 
the emission from ULIRGs, the bulk of the CO luminosity comes from a diffuse inter-cloud 
medium, rather than from the dense gas gravitationally bound in clouds where the stars are 
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formed (Downes et al. 1993; Solomon et al. 1997). The dense gas is better traced by HCN, 
and as shown by Gao et al. (1999) the L FIR / Lh CN -t atio is the same for GMCs to ULIRGs, 
indicating that anywhere in the Universe only the dense gas is relevant to star-formation. In 
ULIRGs the dense/diffuse gas mass ratio (roughly quantified by Lhcn/^co) is particularly 
high which is not a surprise given their merger status and the way the gas responds to a 
merger. 



4.2. Excitation conditions of the molecular gas in J164502 

The velocity /area- averaged brightness temperature ratio between CO J + 1 — > J and 
CO 1 -> is 

Tb(J + 1 — > J) L'(J + 1 — > J) 
TJ+1 ' J T b (l->0) L'(l^O) ' lJ 

Andreani et al. (2000) estimated the CO J = 2 — 1 and J = 5 — 4 luminosities to be 4 x 
10 10 and 7 x 10 9 Kkms" 1 pc 2 , respectively, where we have converted the luminosities to 
the cosmology adopted here. The line ratios estimated from eq. 8 are r 2 i = 0.6 ± 0.2 
and r 54 = 0.10 ± 0.05. The majority of high-redshift systems detected in CO to date are 
strongly gravitationally lensed. This makes estimating the excitation conditions in such 
systems complicated since differential magnification of the high-J lines compared to lower 
transitions may bias the line ratios significantly. There is nothing to suggest that J164502 
is gravitationally lensed, and the above derived line ratios ought to represent the intrinsic 
excitation conditions of the gas, albeit still averaged over the entire galaxy. 

We used the observed line ratios of r 2 i = 0.6 ±0.2 and r 54 = 0.10 ±0.05 as constraints to 
a Large Velocity Gradient (LVG) code in order to gain some insight on the bulk properties of 
the molecular gas in J164502. A lower limit on gas temperature T kin = 40 K is assumed since 
this value is deduced for the bulk of the dust (Dey et al. 1999) and T kin > T dus t is expected 
for FUV-heated gas/dust (e.g. Tielens & Hollenbach 1999). A wide range of temperatures 
(Tkin = 60 — 90 K) offers an acceptable fit to the aforementioned line ratios, with T kin = 70 K 
being the optimal value. This wide range is due to the poor constraints offered by only two 
line ratios while the high gas temperatures may reflect the UV-intense environment of the 
molecular gas in J 164502 and/or turbulent motions heating only the gas, a situation that 
has already been noted for molecular gas in the Galactic center (Rodriguez- Fernandez et al. 
2001). 

A common feature of all the LVG solutions is the low gas densities n(H_2) ~ 300 cm -3 . 
This has already been noted when only the CO (5 — 4)/(2 — 1) line ratio was available 
(Papadopoulos & Ivison 2002) and is a property often found for the bulk of the gas in 
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extreme starbursts (e.g. Aalto et al. 1995; Downes & Solomon 1998). Interestingly for such 
low gas densities all CO J + 1 — > J transitions with J + 1 > 5 have flux density ratios with 
respect to the lowest 1 - transition of S(J + 1 - J)/S{1 - 0) = (J + l) 2 r J+1 j < 0.25 for 
any plausible gas temperature (T^m ~ 40 — 100 K), thus no longer offering the advantage of 
a higher flux density with respect to the latter. 

As a consequence, any estimate of the H 2 mass based on the assumption of optically 
thick, thermalised high-J CO lines could be severely underestimated. Using the 5-4 line 
in J 164502 only, and assuming it is optically thick and thermalised (r 54 ~ 1), one would 
find a molecular gas mass which is an order of magnitude lower than that inferred from 
the J = 1 — line. There is no significant flux advantage between the 2-1 and 1-0 lines, 
since S co (2 - l)/S co {l - 0) = 4T b (2 - l)/T b {l - 0) > 1 even for a very subthermal ratio of 
r 2 i = 0.3, say. This is consistent with the fact that we find the molecular masses derived from 
the two lines to be similar (section 4.1). The above demonstrates the importance of using 
low-J CO lines in order to infer the amount of molecular gas in a galaxy. This is particularly 
true at redshifts beyond 3, since at those redshifts the current mm-interfero meters can only 
hope to detect CO J + 1 — > J, J > 2, and the the excitation bias of the high-J lines could 
become very severe. 



4.3. Comparison with high-redshift sub-mm galaxies 

Although J164502 was originally selected as an ERO and not a (sub) millimetre galaxy, 
subsequent observations at optical/IR and sub-mm wavelengths have shown that it can in 
fact be considered a typical sub-mm galaxy (Dey et al. 1999; see also Smail et al. 1999). 
Our detection of CO J = 1 — in J164502 is therefore the first detection of this transition 
in a galaxy which is thought to be similar to the SCUBA population of dust-enshrouded 
galaxies at high redshifts. J164502 has a far-IR and CO luminosity about ten times that of 
the average values found in ULIRGs, and the amount of molecular gas present in J164502 is 
comparable to the median gas mass (~2x 10 10 M ) of the five SCUBA galaxies detected 
in CO to date (Neri et al. 2003), although a thorough characterisation of the gas content of 
this population has to await CO observations of a large, unbiased sample of submm-selected 
galaxies. Hence, the observations presented here further supports that J164502 is more 
similar in its properties to sub-mm galaxies than a high-redshift analogue of Arp 220. 

Recently, a subject of some debate has been whether the CO emission detected from 
high-redshift sub-mm galaxies originates from a huge, massive resevoir of molecular gas 
~ lOkpc in size or from a much more compact circumnuclear disk, typically of radius r ~ 
100 pc, as seen in local ULIRGs such as Arp 220 (Ivison et al. 2001; Genzel et al. 2003; 
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Downes & Solomon 2003). If, as the observations seem to suggest, the starburst in J164502 
is extended over several tens of kpc that would set apart from the local ULIRG population 
and point toward the former scenario. 
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Fig. 1. — Left: Naturally weighted CO J = 1 — contour map of J 164502 overlaid on the 
UKIRT UFTI K-band image. The resolution of the CO map is 3.1" x 2.8" at PA = 157°, see 
the insert. The contours are -2,2,3,4,5,6, and 7a where a = 0.1 mJy beam -1 . Right: The 
real (top panel) and imaginary (bottom panel) part of the complex visibilities vs. baseline. 
The visibilities have been binned. The solid line represents a point source with flux density 
of l.OmJy. The dashed curves represent a least-square fit of a Gaussian source model to the 
visibilities. 
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Fig. 2.— Left: SCUBA 850-//m jiggle-maps of J164502. The resolution is FWHM = 14", 
and the rms noise level is 1.6 mJy beam -1 . Contours start at 2a. J164502-SMM1 denotes 
a statistically significant detection of 850-/zm emission about 50" south of J164502. Weak 
emission is also seen 20" south-east of J164502. Right: Radial profiles of J164502 and 
J164502-SMM1 (filled circles) compared with the radial profile of the PSF (solid curve). 
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Table 1. VLA Observations 



Date 


Configuration 


Frequency, (GHz) 


Integration time, (hrs) 


2001 Oct 9 


D 


47.235,47.285 


5.2 


2002 Oct 26 


C 


47.235,47.285 


2.9 


2002 Nov 14 


C 


47.235,47.285 


2.4 


2002 Dec 24 


c 


47.235,47.285 


2.6 


2003 Feb 21 


D 


47.235,47.285 


2.6 


2003 Mar 1 


D 


47.235,47.285 


3.8 
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Table 2. Observed properties of J164502 



Parameter 


J164502 


J164502-SMM1 


a (J2000) 


16:45:02.26 


16:45:01.20 


5 (J2000) 


+46:26:26.50 


+46:25:39.0 


ZCO 


1.439 




Jav s co(i-o)<ZV [Jykms- 1 ] 


0.6 ±0.1 


< 0.2 


5850Mm [mjy] 


8 ± 2 a 


6±2 b 


L^ o (l-0) [Kkms-ipc- 1 ] 


(7±1) x 10 10 


< 2 x 10 10c 


Lfir [Lq ] 


9 x 10 12 


7 x 10 12c 



a The average of 7.5/7.9/7.3/9.8 mjy which was the flux densities 
obtained using the four different methods described in section 3. 

b The average of 6.7/6.5/6.0/6.5 mjy which was the flux densities 
obtained using the four different methods described in section 3. 

c This value assumes that J164502-SMM1 has a rcdshift of z = 
1.439. 



